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Ph 1 abstract 

^ I Giant X-ray cavities lie in some active galactic nuclei (AGNs) locating in 

central galaxies of clusters, which are estimated to have stored 10 55 ~ 10 62 erg 
of energy. Most of these cavities are thought to be inflated by jets of AGNs 
on a timescale of > 10 7 years. The jets can be either powered by rotating black 
holes or the accretion disks surrounding black holes, or both. The observations of 
giant X-ray cavities can therefore be used to constrain jet formation mechanisms. 
In this work, we choose the most energetic cavity, MS 0735+7421, with stored 
energy ~ 10 62 erg, to constrain the jet formation mechanisms and the evolution 
of the central massive black hole in this source. The bolometric luminosity of the 
7—i \ AGN in this cavity is ~ 10~ 5 LEdd, however, the mean power of the jet required 

^ | to inflate the cavity is estimated as ~ 0.02LEdd, which implies that the source 

has experienced strong outbursts previously. During outbursts, the jet power and 
the mass accretion rate should be significantly higher than its present values. We 
| construct an accretion disk model, in which the angular momentum and energy 

carried away by jets is properly included, to calculate the spin and mass evolution 
of the massive black hole. In our calculations, different jet formation mechanisms 
£> ; are employed, and we find that the jets generated with the Blandford-Znajek (BZ) 

mechanism are unable to produce the giant cavity with ~ 10 62 erg in this source. 
Only the jets accelerated with the combination of the Blandford-Payne (BP) 
and BZ mechanisms can successfully inflate such a giant cavity, if the magnetic 
pressure is close to equipartition with the total (radiation+gas) pressure of the 
accretion disk. For dynamo generated magnetic field in the disk, such an energetic 
giant cavity can be inflated by the magnetically driven jets only if the initial 
black hole spin parameter a > 0.95. Our calculations show that the final spin 
parameter a of the black hole is always ~ 0.9 — 0.998 for all the computational 
examples which can provide sufficient energy for the cavity of MS 0735+7421. 
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INTRODUCTION 



Giant X-ray cavities have been disc overed in some central g alaxies of cl usters by XMM 
Newton and Chandra [e.g .. Hydra A (|McNamara et al.l 120001 ): RBS 797 ( ISchindler et al. 



20051 )]. which contain huge energy up to 10 



55 



2001|) and MS 0735+7421 flMcNamara et al. 
10 62 erg. They can effectively suppress the cooling of the intracluster medium (ICM) and 
provides a direct evidence for the presence of AGN feedback in galaxy formation and evo - 
lution fe.g. lBirzan et al.l 120041 llVlcNamara et al.ll2005HAUen et al.lbood : ICroton et al.lboOO ). 
It is believed that the cavities are i nflated by the jets launched from the accretion disks in 



active galactic nucleus (AGNs) (see iMcNamara fc Nulsenl 120071 . for a review and the refer 



neces therein). Therefore, the X-ray cavities provide a direct measurement of the mechanical 
energy released by the jets through the work done on the hot gas surrounding them. Mea- 
surements of this energy, combined with measurements of the timescale require d to inflate 



the cavities, can be used to estimate the mean jet power (e.g., lAUen et al.ll2006l ) 



The jet formation mechanisms in AGNs have been extensively studied in the last several 
decades. The currently most favored jet formation mechanisms contain two categories, i.e., 



the Blandford -Znajek (BZ) process (jBlandford &: Znajekl 119771 ) and the Blandford-Payne 
(BP) process ( Blandford fc Payne 19821 ). In the BZ process, energy and angular momen- 



tum are extracted from a rotating black hole and transferred to a remote astrophysical 
load by open magnetic field lines. In the BP process, the magnetic field threading the 
disk extracts energy from the r otating gas in the accretion disk to power the jet/outflow. 



The hybrid model proposed by iMeierl (119991 ). as a variant of the BZ model, which com- 



bines the BZ and BP effects through the large scale magnetic field threading the accre- 
tion disk outside the ergosphere and the rotating plasma within the ergosphere. The 
ordered large-scale magnetic field threading the disk is a crucial ingredient in jet forma- 



scale with the gas/radiation pressure of the accretion disk 


(e.g.. Moderski & Sikora 1996; 


Ghosh & Abramowicz 


1997: Livio et al.l 


1999: 


Armitaee 


k Nataraian 1999: 


Nemmen et al. 


2007: Wu & Cao 2008: 


McNamara et al. 


2009: 


Wu et a . 


201 


.1). Livio et al. 


(1999) pointed 


out that even the calculations of Ghosh & Abramowicz ( '. 


L997 


) have overestimated the power 



of the Blandford-Znajek process, since they have overestimated the strength of the large-scale 
field threading the inner region of an accretion disk, and then the power of the Blandford- 
Znajek process. The strength of the large-scale field scales with the disk thickness, and it 
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is very weak if the field is created by dynamo processes for thin disk cases (ITout fc Pringle 



19961 ; iLivio et al.l Il999l ). The maximal jet power extracted from an accretion disk can be 



estimated on the assumption that the toroidal field c omponent is of the same order of th e 
poloidal field component at the disk surface (see, e.g.. lLivio fc Pringldll992l ; lBiskamplll993l ). 



The giant cavity in MS 0735+7421, one of the most energetic outbursts in AGNs, which 
stores a total ~ 10 62 erg of ene rgy with the timescale required to inflate the cavity ~ 10 8 years 

10 46 erg s" 1 



(McNamara et al 



(IMcNamara et al 



2005 



2005, 2009). It implies that the mean jet power should be 



Gitti et al.l 120071 ). which is about two percent of the Eddington 



l uminosity (L FlM ) of a 5 x 1O 9 M (M is the solar mass) black hole in MS 0735+7421 
(IMcNamara et al.ll2009l ). However, the central AGN in this cavity shows very low optical 
nuclear emission (Li < 2.5 x 10 42 erg s _1 (Li is the luminosity in I band), which implies its 
accr etion disk accreting at the curr ent rate being too low to power strong jets with ~ 10 46 erg 



McNamara et al.l (120091 . l201ll ) argued that the BP process solely is insufficient to power 



the giant outbursts in this source. They suggested that the jet is alternatively powered by 
the rotational energy of a extremely Kerr black hole (~ 10 62 erg ) , whic h roughly corresponds 



to the total energy stored in this cavity. However, ILivio et al.l (Il999l ) pointed out that the 
magnetic field strength near the horizon of a rotating black hole is dominantly determined 
by the structure of the inner region of the disk. Thus, it is still quite doubtful if such a 
faint accretion disk can maintain a strong field near the horizon to power powerful jets with 
~ 10 46 erg s- 1 . 

It is well known that quasars are powered by accretion onto massive black holes, and the 
growth of massive black holes could be dominantly go verned by m ass accretion in quasars. 
The massive black holes are therefore the AGN relics (ISoltanlll982l ). The faint AGN in MS 
0735+7421 contains a massive black hole with mass of 5 x 10 9 M Q , which may probably have 
experienced a quasar phase before. During its quasar phase, the black hole was accreting at 
relatively high rates, and therefore powerful jets can be magnetically accelerated from the 
region near the black hole either by the BP or BZ precesses, or both. In the same period of 
time, the central massive black hole grew up, and the hole was spun up simultaneously. In 
this work, we choose the most powerful cavity, MS 0735+7421, to constrain the jet formation 
mechanisms and the mass/spin evolution of the massive black hole in this source. 
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MODEL 



2.1. The equations of accretion disk 

The faint AGN in MS 0735+7421 may probably have experienced a quasar phase before. 
During its quasar phase, the black hole was accreting at relatively high rates, and therefore 
its ac cretion disk can be described by the standard thin disk model (jShakura fc Sunyaev 
19731) . In this work , we consider a relatiy istic thin accretion disk around a Kerr black hole 
(jNovikov fc Thorne 1973 ; Manmoto 2000 ). and the metric around the black hole reads (the 
geometrical unit G = c = 1 is adopted): 



t 2 /S. A t 2 

dt 2 + —(d(f)- udt) 2 + —dr 2 + dz 2 



A r 2V r A 

A = r 2 - 2Mr + a 2 , 
A = r 4 + r 2 a 2 + 2Mra 2 , 



UJ 



2Mar 

j 

where M is the mass of the black hole, J and a are the angular momentum and specific 
angular momentum of the black hole respectively, and u is the dragging angular velocity of 
the metric. 

The model of the geometrically thin, opti cally thick and Kep l erian accretion disk sur- 
rounding a Kerr black hole was developed by iNovikov fc Thornd (119731 ). The structure of 
the accretion disk may be altered in the presence of the large-scale magnetic field thread- 
ing the disk, which can drive the outflows/jets from the disk, and a fractio n of the angular 
mom e ntum and energy of the disk is carried away in the outflows/jets (e.g., lOgilvie fc Livio 



19981 ; ICao fc Spruitl |2002| ). In order to explore the magnetically accelerated jets from the 



accretion diks, we properly consider the effects of the jets on the accretion disk. The basic 
equations of a relativistic thin disk with magnetically driven outflows/jets are summarized 
as follows. 



The continuity equation is 



M = -2ttA 1/2 E 



(2) 



where v r is the radial velocity of the accretion flow, S = 2pH is the surface density , and M 
is the mass accretion rate. As we consider fast moving outflows/jets in this work, the mass 
loss rate in the ouflow/jet is neglected in the continuity equation of the accretion flow. 
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The radial momentum equation is 

O+VO - (). \ 

+ 9m = 0, (3) 



7^M (fl-fl+)(Q-Q- 



r 4 A fi+ft- 



where Q is the angular velocity, and the Lorentz factor 7^ of the rotational velocity v$ is 
given by 



14> 



and Q = Q — u. In this work, we consider the geometrically thin disks, and the gas pressure 
gradient in the radial direc tion can be omitted in the radial momentum equation of the disk 
(jShakura &: Sunyaevlll973l ). The first term in Equation represents the net force exerted 
on the gas moving with a circular angular velocity Q in the Kerr metric. It should be balanced 
with the radial magnetic force g m induced by the large-scale magnetic field threading the 
accretion disk without considering the radial pressure gradient in the thin disk case. The 
Keplerian angular velocities of the prograde (+) and retrograde (— ) motions are 

nt 



k r 3 / 2 ± aM 1 / 2 ' 
and the radial magnetic force is 

g m = B t B z /2ttE, (4) 

where B r and B z are the radial and vertical components of the magnetic fields at the disk 
surface, respectively. The inclination of the field line at the disk surface k is defined as 
^0 = B z /B r . 

The angular momentum equation is 

_|£^ + ^ W)+ T m r = 0, (5) 
2/K ar ar 

where the angular momentum of the accretion flow L is 

L — ) 

r 

and the height-integrated viscous tensor is 

A 3/2 A l/2 3 

w; = a -g —W, 

where a is the Shakura-Sunyaev viscosity parameter. The height-integrated pressure W = 
2HP tot , where the total pressure -Ptot — Pgas -Frad Pgasi -frad ^rid P ai — (1 — 0)Ptot &re 
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the gas pressure, radiation pressure and magnetic pressure respectively (1 — (5 is the ratio 
of the magnetic pressure to the total pressure). The scale height H of the accretion disk is 
given by 

H 2 = C y/(L 2 -a% (6) 



where c s = y-Ptot/p is the sound speed of the gas in the disk. 

The first two terms in Equation (jSJ) represent the change rate of angular momentum 
and the transfer rate of angular momentum caused by the viscous tensor in the gas, and 
the third term is the rate of angular momentum carried away by the magnetically driven 
outflows/jets. The magnetic torque exerted on the accretion flow due to the outflows/jets is 

_ B Z B V R 
Ztt 

where B v is the toroidal component of the field strength at the disk surface. Equation (J7J) is 
the vertical integrat ion of magnetic torque in the ideal MHD angular momentum equation 



f lCao fc Spruitl 120021 ) . 



The energy equation is 



r 6 \dr J 3kS 



where v is the viscosity coef ficient, uT^dfl/dr = —aW /r in a- viscosity, and T is the temper- 



ature of the gas in the disk (lAbramowicz et al.lll996l ). The left term in this equation is the 



surface heat generation rate caused by turbulence of the gas in the disk, which is balanced 
with the cooling rate of the disk. The radial advection of energy is neglected in the thin disk 
case. The opacity R of the gas is given by 

k = K es + n s = 0.4 + 0.64 x 10 23 pT~ 7/2 cm 2 g~\ 
where K es and Kg are the electron scattering opacity and free-free opacity respectively. 



2.2. The jet formation mechanisms 

In the general form of the BZ mechanism, the jet power Lbz can be estimated with 
(IGhosh fc Abramowiczlll997l ) 

Lbz = ^oo 2 F Blr 2 u (J/J max ) 2 c, (9) 

where rn is the horizon radius, B± is the component of the magnetic field normal to the 
black hole horizon, J and J max = GM 2 /c are the angular momentum and maximum angular 
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momentum of a black hole, and Up = fip^H — ^f)/^h * s a factor at the black hole horizon 
determined by the angular velocity of black hole and that of the magnetic filed lines. In this 
work, we simply adopt B ~ B (B 2 = B 2 +B 2 „), which is the maximal magnetic field strength 
in the disk (see lGhosh & Abramowiczl 119971 . for the details). It is easy to conclude that the 
maximal jet power L^% x corresponds to = 1 /2Q-R. The power of the jets accelerated from 



an accretion disk can be calculated with (e.g., iLivio et al.lll999t ICaol 120021 ) 



rout g £ 

p rQ2nrdr, 



4tt 



(10) 



where B p (B 2 = B 2 + B 2 ) is the poloidal field strength, and r in and r c 
outer radius of the accretion disk respectively. 



are the inner and 



Tout fc Pringld ( 119961 ) suggested that the typical size of the magnetic fields produced 
by dynamo processes is roughly around the disk thickness H. The large-scale field can 
be produced from the small-scale field created by dynamo processes as B(X) oc A -1 for 
the idealized case, where A is the length-scale of the field. For the magnetically launching 
problem for jets, the size of magnetic field lines > R is required in order to accelerate the 
jets/outflows efficiently. Thus, ILivio et al.l (Il999l ) proposed that the poloidal magnetic fields 
B should be ~ (H/r)B if they are generated through dynamo processes in the accretion 
disk. The jet power of the BZ and BP mechanisms is given by 



1 / H \ 

Lbz = ( B ±— ) r n( J / 4ax) 2 c, 

61 V r / max 



and 



r ° ut B B ( H 

Ait \ r 



dr, 



(12) 



on the assumption of magnetic field generated through dynamo processes, where (B±H/r) max 
is the maximal (B±H/r) in the inner region of the accretion disk. 

The difficulties in estimating the jet power, either for the BP or the BZ processes, are the 
magnetic field strength and the field geometry, which are still highly uncertain. The conven- 
tional estimates of the field strength are more or less based on the assu mption of equiparti- 
tion between magnetic pressure and gas pressure /radiation pressure fe.g..lModerski fc Sikora 



1996 



2007 



Ghosh fc Abramowicz 



Wu k Cao 



2008 



1997; 



Livio et alJl999 



McNamara et al. 2009; 



Armitage fc Natarajanlll999i ; lNemmen et al 



Wu et al. boil ). In this work, we adopt a pa- 
rameter (3 to describe the field strength in the disk. The magnetic pressure P m = (1 — (3)P tot 
(-Ftot = -Pgas + -Frad); an d therefore the field strength B = a/(1 — /3)8nP tot . In principle, the 
magnetic field strength and the structure of the accretion disk are available by solving a set 
of dynamical equations of the disk described in Section 2.1 simultaneously, if the value of (3 
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and the field geometry are known. Howev er, the geom etry of the magnetic field of the disk 
is in principle a global problem (see, e.g., ISpruitl l201ll for the detailed discussion, and the 
references therein). It depends not only on the ini tial strength and geometry of the magnetic 



field as indicated by some M HD simulations (e.g.. Ilgumenshchev et al.ll2003l ; iBeckwith et al. 



20081: iMcKinnev et al. 



disks (ILubow et al.lll994 : 



20 



2j), but also on the processes of d iffusion and advection in accretion 



Cadl201ll ; lMcKinney et al.ll2012l ). which is beyond the scope of this 



paper. The inclination of the field line at the disk surface plays a key role in launching jets. 
More specifically, the angle of field lines inclined to the mid-pla ne of the disk is required t o 
be less than 60° for launching jets from a Keplerian cold disk (IBlandford &: Paynd Il982l ). 
This critical angle can be slightly larg er than 60° if the internal energy of the gas in the 
disk is considered ( ICao fc Spruitl I1994J ) . This critical angle could be larger than 60° even 
for the cold gas magnetically lau nched from the inner edge of the accretion disk very close 
to a rapidly spinning black hole (ICaolll997l ). which indicates t hat the spin of the black hole 
may help to launch jets centrifugally by cold magnetized disks ( ]Caolll997t ISadowski fc Sikora 
20101 ). This is because the frame drag effect of the spinning black hole may help to accelerate 
the outflows. For realistic cases, the internal energy of the gas is included in the Bernoulli 
equation of the outflow, and an additional force due to the pressure gradie nt helps to ac- 
celerate the gas in the outflow (e.g., see equation 14 in ICao &: Spruitl 1 19941 ). In this case, 
the outflow can be launched if the inclination angle is slightly larger than the critical value. 
As the final derived accretion disk structure is insensitive to the precise value of the field 
inclination at the disk surface, we simply adopt the inclination of 60° in all our calculations. 
The strength of the azimuthal component of the field at the disk surface is mainly deter- 
mined by the properties of the magnetically jets. The model properly including both the 
accretion disk solution and jet solution is very complicated, which is beyond the scope of 
this work. It was pointed out that the fast moving jets/outflows always correspond to the 
case B v <C B p , which B ip ~ B p is satisfied for slowly moving outflows with relatively high 



mass loss rate (ICao fc Spruitl 1 19941 ; lOgilvie fc Liviolll998l : ICao & Spruitll2002l : ICaoll2002h . In 

the case of the problem considered in this work, we focus on the fast moving jets from the 
disk, and B v <C B p should be satisfied. In this work, we adopt £^ = 0.1 (B v = £ipB p ) in 
most of our calculations. 



In summary, we calculate the jet power L\ e t with two models in this work, i.e. 



Model A (general BZ+BP mechanisms): Lj et = Lbz + Lbp, where Lbz and Lbp are 
calculated with Equations and (ITU]) respectively; 



Model B (Livio's model): Lj et = Lbz + Lbp, where Lbz and Lbp are calculated with 



Equations ( ITT]) and ffl~2"]) respectively. 



The main difference between model A and B is the estimate of the field strength in the 
accretion disk. Finally, the energy released through BZ+BP mechanism during the outbursts 
is available with t 

E to t = E BZ + E BP = / Lj Ct dt, (13) 
Jo 

where t is the duration of the AGN outbursts. 



2.3. The evolution of the black hole 



For a massive black hole s urrounded by an accretion disk, the blac k hole mass and spin 
evolution is described by (e.g., iModerski fc Sikoralll996l ; iLu et al.lll996l ) 



dt M Ums ms) Mc 2 \kh h 



2a 



(14) 



dlnM M 



Lbz 



dt M~ LL1 ~ (1:J) 
where a (0 < a < 1) is the dimensionless angular momentum of black hole, M = dM/dt is 
the accretion rate, M is the mass of black hole, j ms and e ms are the dimensionless specific 
angular momentum and energy of the gas in the accretion disk at the marginally stable orbit 
respectively, i\ is the dimensionless angular velocity at the horizon of the black hole, Lbz 
is the jet power due to the BZ process, and k = flp/Qy, < 1 is a constant, where f2p is the 
angular velocity of magnetic field lines at the horizon. 

The time evolution of mass accretion rate is still quite unclear, which may be dependent 
on the circum-nucl ear gas near the black hole or /and feedback of the quasar at the center 
of the galaxy (e.g., |Pi Matteo et al.l 120051 ) . The observed Eddington ratio distribution for 
AGNs can be roughly described by a power-law distribution, which is consistent with the 
self-regulated black hole growth model. In this model, the feedback of AGNs produces a 
self-regulating decay or blowout phase after the AGN reaches som e peak luminosity and 
begin s to expel gas and shut d own a ccretion ( Hopkins et al. 2005al lbl; lliopkins fc Hernquist 



20091 ). iHopkins fe Hernquistl (120091 ) suggested the quasar luminosity (then the accretion 
rate) evolving with time can be described by 
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where m = M/M-^dd is the accretion rate in units of Eddington accretion rate (Msdd = 
1.5 x 10 18 M/M Q g s -1 is the Eddington accretion rate), tq is a consta nt and m n is the peak 



accret ion rate of quasar. The index f3 L ~ 0.6 is adopted as suggested by lHopkins fc Hernquist 
(120091 ). We define the lifetime tq of quasar as 

o.oi / ^ \ -fa 



tc 



dt = - 



m 
rh Q 



exp 



) dlogm, 

mo J 



(17) 



where the integr al upper limit m 
line AGNs (e.g. 



0. 01 is adopted as t he typical lower limit on broad- 
Kollmeier et al.l 120061 ; iTrump et al.l 120111 ) . We assume m = 1 in all our 
calculations. The quasar lifetime £q ~ 10tq for /3l = 0.6. 



3. RESULTS 

As described in Sections 2.1 and 2.2, the structure of an accretion disk surrounding 
a rotating black hole with mass M and spin parameter a can be calculated by solving 
Equations (J2J) — (JS]) for different jet formation mechanisms, provided the values of the disk 
parameters, a, (3, and M, are supplied. Based on the derived disk structure and jet power, 
the mass and spin evolution of the black hole can be calculate with Equations (fT4"|) and ffT5"]) . 
if suitable initial conditions and the mass accretion rate as a function of time are specified. 
The conventional value of viscosity parameter a = 0.1 is adopted in all our calculations. In 
Figure dj we find that the disk structure has been significantly altered by the magnetically 
accelerated jets. The temperature of the accretion disk with outflows/jets is obviously lower 
than that of the standard thin accretion disk, because a fraction of the gravitational energy 
released in the disk is carried away in the outflows/jets. The temperature of accretion disk 
decreases significantly for model A than that for model B, because the jet power is higher 
for model A if all other model parameters are fixed (see Figured]). 

The central AGN in the cavity of MS 0735+7421 contains a 5 x 1O 9 M black hole, 
which shows very low optical nuclear emission (Lj < 2.5 x 10 42 erg s -1 in I band). The 
model calculations for this source should satisfy three observational quantities, total energy 
stored in the cavity ~ 10 62 erg, present mass of the cent ral black hole 5 x 10 9 M C ? ) , and the 



timescale of the jet activity in this AGN ~ 10 8 years (IMcNamara et al.l 12005k 120091 ). It 
implies that the giant cavity in this source was inflated by the jets when the black hole was 
accreting at relatively high rates previously. In all our calculations, the initial black hole mass 
Mq is chosen in such a way to let the final black hole mass be 5 x 1O 9 M as that in the AGN 
of MS 0735+7421. Thus, the time-dependent mass accretion rate m{t) is a crucial ingredient 
in our model calculations. In this work, we adopt a power- law time dependent m (t) with an 



exponential cutoff at the high end rho suggested by lHopkins fc Hernquistl (120091 ). It is found 
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that the mass accretion rate declines very rapidly for this time-dependent rh(t) when tq is 
small, i.e., a short quasar lifetime (see Figure [2]). 

The total energy released during the outbursts as functions of time with different values 
of the model parameters for different jet formation models are plotted in Figures [3J171 In 
Figures [3] and HI we plot the energy released during the outbursts as functions of the outburst 
time with different sets of the model parameters for different initial black hole spin parameters 
a = 0.1 and 0.95, respectively. It is found that the total energy released in the cavity is 
higher for model A than that for model B. The total energy increases with tq, and much 
energy is carried away in the jets provided the initial black hole spin parameter is high, if 
all other parameters are fixed. 

The black hole mass and bolometric luminosity of broad-line AGNs can be well estimated 
with the observations of the SEDs and broad-line width, and therefore the Eddington ratio 



are around 0.1 — 0.4 (e.g. 


McLure & DunloD 


Kauffmann & Heckman 


2009; 


Shen k Kellv 



20121 ). The time evolution of accretion rate is still quite uncertain. A constant accretion rate 



is usually assumed for AG Ns in most of the previous studies of the cosmological evolutio n 



of massive black holes (e.g., lYu &: Tremaineil2002l ; iMarconi et al.ll2004l ; IShankar et al.l 120041 ) . 



The jet power is dependent on the mass accretion rate rh, and for comparison, we alterna- 
tively adopt a constant mass accretion rate, fn = 0.2 5, which is around the peak value of the 
distribution in the sample of iKollmeier et al.l (120061 ). to calculate the total energy released 
in the cavity as functions of time (see Figure [5]) . 

We compare the relative importanc e of the BP and BZ mechanisms in Figure El Similar 
to that suggested by iLivio et al.l (119991 ). we find that the BP power always dominates over 
the BZ power even if the black hole is rotating rapidly. The results calculated with different 
values of /3 are plotted in Figure [71 The jet power decreases significantly with increasing 
(3. In most of our calculations, (3 = 0.5 is adopted, which is almost the upper limit on the 
magnetic field strength, and then the jet power. 



The evolution of black hole spin for different jet formation models is given in Figure [HI 



4. DISCUSSION 

The physics of jet formation is still quite unclear. As described in the previous sections, 
the total energy released by the jets of the AGN and the mass/spin evolution of the central 
massive black hole can be calculated with suitable initial conditions for a given jet formation 
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model. The most powerful outbursts observed in MS 0735+7421 provide constraints on 
different jet formation models. 

The basic physic s of the jet formation model s ado pted in thi s work is the same as 
the previous works of I Ghosh fc Abramowica (119971 ) and iLivio et al.l (119991 ) , however, their 
estimates of the magnetic field strength of the disk are based on the conventional accretion 
disk models without magnetic field. This is a good approximation in the weak magnetic 
field the disk structure has not been altered significantly by the field, while the 

assumption becomes invalid for strong field cases. In this work, we calculate the structure of 
the accretion disk including the magnetic torque exerted on the disk and the energy carried 
away by the jets (see Section 2.1). 

The giant cavity in MS 0735+7421 stores a total ~ 10 62 erg of energy, which is the 
most energetic outbursts of AGN so far discovered, if the cavity is inflated by the jets of the 
AGN. The observations of MS 0735+7421 provide very useful information that can constrain 
the detailed jet formation models. A power-law time dependent m(t) with an exponential 
cutoff at the high end is adopted in most of our calculations, which is derived based on 
the scenario of the feedback of AGNs produ cing a self- regulating de cay or blowout phase 



after the AGN reaches some peak luminosity (iHopkins et al.ll2005al Jbl; iHopkins fc Hernquist 



20091 ). It is found that the total energy released in the jets can be as high as > 10 62 erg only 
in the calculations with model A, i.e, the strength of large-scale magnetic field fields scales 
directly with the total pressure of the accretion disk, if the initial black hole spin parameter 
a = 0.1 is adopted (see Figure [3]). The another jet formation model seems not to be able to 
produce such a giant cavity containing ~ 10 62 erg even if equipartition between magnetic and 
total(radiation+gas) pressure in the disk is assumed (i.e., (3 = 0.5). As discussed in Section 
2.2, <C B p is required for the fast moving jets from the disk, and we adopt £ v = 0.1 
(Bip = iipBp) in most of our calculations. For comparison, we also plot the results calculated 
with B v = B p in Figures [3] and HI We find that the maximal jet power can be extracted with 
model A is insensitive to the value of £ v adopted, because most of the gravitational energy 
released in the accretion disk is transported into the jets by the magnetic field. The results 
of the calculations with £^ = 0.1 are given in all other plots. 

The initial black hole spin parameter a can be higher than 0.1, for example, a rapidly 
rotating black hole can be formed after a mer ger of the black hole binary with mass ratio 
approaching unity ( {Hughes fc Blandfordl 120031 ) . Thus, we also calculate the problem with a 
relatively high initial black hole spin parameter = 0.95 in Figure HI It is not surprising 
that the observed feature of the cavity can be reproduced by the calculations with model 
A, as much energy can be extracted for a faster rotating black hole. One can find that 
the total energy released in the jets can be as high as ~ 10 62 erg if qq = 0.95 and tq = 
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10 8 years, even for model B, i.e., the energy density of the large-scale magnetic field is 
scaled with p tot (H / r)~ 2 (see Equations [11] and [12] in Section 2.2), based on the scenario 
that the strength of the magnetic field generated by dynamo processes with typical size A 
decays with A -1 (ITout fc Pringld Il996l ) . If this is the case, the life timescale of this AGN 
tq ~ IOtq = 10 9 years, which is roughly consistent with the estimates in the previou s 
works (e.g., bi Matteo et al.ll2005[ iHopkins k Hernquistl[2009t kellv et alJboiot Icaolboioh . 
This implies that the duty cycle of radio activity is around 0.1, if model B is responsible 
for the jet formation in this source. Besides the time- dependent mass accretion rate rh(t), 
we alternatively adopt a constant accretio n rate in our calcula tions, rh = 0.25, which is 
the typical value of broad- line AGNs (e.g., iKollmeier et al.l 120061 ) . We find that the results 



change little if a different value of mass accretion rate in the range of 0.1 — 0.5 is adopted. 
The results obtained with a constant accretion rate are similar to those carried out with the 
time-dependent mass accretion rate, i.e., the jets can provide sufficient energy to inflate the 
cavity of MS 0735+7421 within the timescale of 10 8 years with model A or model B if an 
initial black hole spinning at ao = 0.95 (see Figure [5]). 

We compare the relative importanc e of the BP and BZ mechanisms in Figure Similar 
to that suggested by iLivio et al.l ( 119991 ). we find that the BP power always dominates over 
the BZ power even if the black hole is rotating rapidly. The reason is that the area of the 
disk launching the jets is significantly larger than the surface area of the black hole, and the 
strength of the field th reading the horiz on of the hole is comparable with that in the inner 
region of the disk (see ILivio et al.lll999l . for the detailed discussion). The initial black hole 
spin ao can change the power output of both the BZ and BP mechanisms as shown in Fig. 
[6] The jet power produced by BZ mechanism can vary over one order of magnitude, which i s 
consistent with the results of MHD simulation for a thin disk in iTchekhovskoy et al.l ( 120101 ). 
Black hole spin can play a m ore important role in the jet formation for thicker accretion disks 
( ITchekhovskoy et al.l 120101 ) . The large- scale magne tic field is dragged in by the advection 
dominated accretion flow (ADAF), and ICaol ( 120111 ) found that the magnetic field strength 
of the flow near the black hole horizon can be more than one order of magnitude higher 
than that in the ADAF at ~ 6GM/c 2 , due to its large radial velocity of the accretion near 
the black hole horizon. It is still unclear if such effect is important in the thin accretion 
disk, which is beyond the scope of this work. We compare the results with different values 
of (3 adopted in Figure The jet power decreases significantly with increasing (3. In most 
of our calculations, (3 = 0.5 is adopted, which is almost the upper limit on the magnetic 
field strength, and then the jet power. Our main conclusions will not be altered even if a 
lower (3 is adopted, while the total energy derived with model B may be insufficient for the 
energy stored in the cavity of MS 0735+7421 provided a lower (3 is adopted. Our results 
show that strong magnetic field nearly equipartition with the total pressure (radiation+gas 
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pressure) in the accretion disk is indeed required at least for the most energetic cavity of MS 
0735+7421, which provides a useful clue on the formation of the large-scale magnetic field 
in the accretion disk. For other less energetic cavities, the magnetic field strength can be 
significantly lower than the equipartition value, and less strict constraints are set on the jet 
formation mechanisms. 

It is found that the final spin parameter Of is always very high (af > 0.9) for all the jets 
that can provide sufficient energy for the cavity of MS 0735+7421 (also see Figure E]), either 
for a low or high initial spin parameter a (af = 0.923 and 0.998 for a = 0.1 and a = 0.9 
respectively). This means that the central black hole of the AGN in MS 0735+7421 should 
be rotating very rapidly now. 
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Fig. 1. — The disk temperature of a relativistic thin accretion disk with magnetically driven 
outflows, where a = 0.95, (3 = 0.5 and m = 0.25 are adopted. The black and red lines are 
for the results calculated with Model A and Model B respectively. We also plot the results 
of a relativistic thin accretion disk without outflows for comparison (green lines). 
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Fig. 2. — The evolution of accretion rate as functions of outburst time. The colored lines 
correspond to different values of tq, tq = 5 x 10 7 (black), 10 8 (red), and 10 7 years(green), 
respectively. In all our calculations, the initial black hole mass M is chosen in such a way 
to let the final black hole mass be 5 x 1O 9 M . 
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Fig. 3. — The energy released during the outbursts as functions of the outburst time with 
do = 0.1 and (3 = 0.5. The different color lines represent the results calculated for different 
jet formation models and different values of azimuthal component of the field at the disk 
surface (see the Figure). The solid, dashed and dotted line are for the results with tq = 10 8 , 
5 x 10 7 , and 10 7 years, respectively. 
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Fig. 4. — The same as Fig. [3j except that the initial spin a 



= 0.95 is adopted. 
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Fig. 5. — The energy released during the outbursts as functions of outburst time with the 
initial spin parameter, a = O.l(dashed) and 0.95(solid), respectively. The model parameter 
(3 = 0.5, £ v = 0.1, and a constant accretion rate m = 0.25 are adopted. The colored lines are 
for the results calculated with different jet formation models, Model A(black) and B(red). 
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Fig. 6. — The energy released by BZ+BP mechanisms (model A) as functions of outburst 
time, where (3 = 0.5, £ v = 0.1, and tq = 1 x 10 7 years are adopted. The black and red lines 
are the results calculated for the BP and BZ mechanisms respectively. The solid and dotted 
lines correspond to different initial conditions, a = 0.95 and a = 0.1, respectively. 
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Fig. 7. — The energy released during the outbursts as functions of outburst time with differ- 
ent values of (3 for model A, where £^ = 0.1, tq = 10 7 years are adopted in the calculations. 
The colored lines are for the results calculated with (3 = 0.99(green), 0.9(black), and 0.5(red), 
respectively. The solid and dotted lines correspond to different initial conditions, ao = 0.95 
and a = 0.1, respectively. 
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Fig. 8. — The evolution of black hole spin as functions of outburst time for model A, where 
f3 = 0.5 and ^ = 0.1 are adopted. The colored lines correspond to different values of tq, 
tq = 5 x 10 7 (black), 10 8 (red), and 10 7 years(green), respectively. The dashed and solid lines 
correspond to different initial conditions, ao = 0.1 and ao = 0.95, respectively. 



